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Abstract 

Critical phenomena in non-equilibrium systems have been studied by means of a 
wide variety of theoretical and experimental approaches. Mode-coupling, renormaliza- 
tion group, complex Lie algebras and diagrammatic techniques are some of the usual the- 
oretical tools. Experimental studies include light and inelastic neutron scattering. X-ray 
photon correlation spectroscopy, microwave interferometry and several other techniques. 
Nevertheless no conclusive treatment has been developed from the basic principles of a 
thermodynamic theory of irreversible processes. We have developed a formalism in 
which we obtain correlation functions as field averages of the associated functions. By 
applying such formalism we attempt to find out if the resulting correlation functions 
will inherit the mathematical properties (integrability, generalized homogeneity, scaling 
laws) of its parent potentials, and we will also use these correlation functions to study the 
behavior of macroscopic systems /ar/rom equilibrium, specially in the neighborhood of 
critical points or dynamic phase transitions. As a working example we will consider the 
mono-critical behavior of a non-equilibrium binary fluid mixture close to its consolute 
point. 

Keywords: dynamic critical phenomena, extended irreversible thermodynamics. Corre- 
lation function. Critical exponents. Binary fluid mixture. 



1 Introduction 

The physical description of systems close to critical points has been attempted from sev- 
eral standpoints ranging from the pure description of the critical opalescence of liquids by 
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Andrews in the XlXth century to modern approaches based on the application of Quantum 
Field theoretical tools such as various kinds of renormalization procedures [fTl|, non-affine 
algebraic structures, etc. Thermodynamically based descriptions of such systems are also 
available [l2l[3l|4l. From the days of van der Waals, passing through Landau and Ginzburg, 
thermodynamic studies of critical properties have been related to geometrical properties of 
the thermodynamic state space. Different geometrical features were invoked: accidental geo- 
metrical relationships to describe coexistence or multi-criticality [3J; a Riemannian geometric 
conjecture relating the thermodynamic curvature with the correlation length and the free en- 
ergy [El, etc. 

Since critical phenomena involves strong dynamic behavior of fluctuations which induces 
abnormal transport, a non-equilibrium study will be adequate. Most dynamic studies of criti- 
cal systems have been carried out based on the development of a description of internal modes 
known as mode-coupling theory. The pioneering works of Fixman [6J and Kawasaki Q as 
well as others |[8l HI [TOl [IB developed into a robust theory compatible with the experimen- 
tal evidence from scattering and sound attenuation studies [[T2l [T3l [T4l [T51 . Unfortunately 
the mode-coupling framework commonly invokes the use of the so-called generalized hydro- 
dynamics [[T6ll based upon a linear (classical) version of the thermodynamics of irreversible 
processes. Linear irreversible thermodynamics assumes that the media under consideration 
could be taken as an equilibrium thermodynamic system pointwise, that is, that the quantities 
used to characterize an equilibrium state exist as continuous space-time functions. It is im- 
portant to stress that not all thermodynamic systems could be viewed as these continua since 
in most cases there are domains or regions which either being composed of aggregates (like 
colloids and macromolecules on a suspension) or because of an uneven distribution of some 
property (anisotropic stress distribution, concentration domains, etc.) violate this pointwise 
equilibrium hypothesis which has been called the local equilibrium hypothesis fil7il . 

As it is known linear irreversible thermodynamics involves an assumption in the form 
of the already mentioned local equilibrium hypothesis stating that the description of sys- 
tems close to equilibrium depends just on a local field evaluation of equilibrium quantities. 
Notwithstanding this condition mode-coupling theory often includes non-local effects due 
to the presence of non-equilibrium phenomena in the form of viscous and diffusive modes 
[[T41 . These effects induce non-locality in space or time because of the presence of regions or 
domains. For example, in a system with concentration domains, mass diffusion occurs and 
so the concentration in one point is changing in time, thus this point is not on an equilibrium 
condition. If the concentration gradients are small, as is the case when the equilibrium state 
is almost reached, one can find regions in which the concentration field is almost homoge- 
neous and, in this case the local equilibrium will be a good approximation. Of course, this 
is not the case if the system presents very large values of fluctuations, since the relaxation 
of these fluctuations will induce domains that will produce transport processes, violating the 
local equilibrium assumption. 
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In order to deal with cases in which the local equilibrium hypothesis fails, several theories 
usually englobed under the name extended irreversible thermodynamics have been developed. 
The common issue in all these extended irreversible thermodynamic theories is the inclusion 
of the non-local effects in the form of additional thermodynamic forces or, in some cases the 
addition of higher order terms in the gradients of the non-conserved properties. For example, 
if there is an uneven concentration distribution, this fact will induce mass transport. Mass 
transport as a non-equilibrium process will affect the whole thermodynamic description of 
the system, so that, for example its thermal diffusivity, its electrical conductivity and its heat 
capacity will change [fTSl . 

Due to the afore mentioned reasons in this work we attempt to introduce an extended irre- 
versible thermodynamical description of the behavior of critical systems based on properties 
of the extended thermodynamic space (integrability, uniform stability, generalized homo- 
geneity, etc.) that takes into account these large fluctuation non-local effects and in this way 
analyze the relation between criticality and dissipative phenomena. 



2 Thermodynamic formalism 



Following the tenets of Extended Irreversible Thermodynamics (FIT) we shall start our dis- 
cussion by assuming that a generalized entropy-like function ^ exists, which may be written 
in the form [,191 : 



udU dv ^ ^dCi ^ d^, 

—;— = T \— — Vv— > Uj—, / '■ —r-\ (1) 

dt ^ dt ^ dt dt ^ ^ dt ^ 

or as a differential form 

i j 

Equation [2]resembles the proposal made by Chen and Fu [|20ll that 

Tdt^ = dtU + pdtV - J2 f^UtCi + dt^j (3) 

i j 

We see that equations ([2] and [3]) are nothing but the formal extension of the celebrated Gibbs 
equation of equilibrium thermodynamics for the case of a multi-component non-equilibrium 
system. The quantities appearing therein are the standard ones,i.e. T is the local temperature, 
p and V the pressure and volume, etc. Xj and $j are extended thermodynamical fluxes 
and forces. These extended forces and fluxes are the new elements of FIT, the ones that 
take into account the aforementioned non-local effects. In the case of a binary fluid mixture 
close to its critical consolute point we will take our set of relevant variables to consist in 
the temperature T(r, t) and concentration of one of the species C2(r, t) as the slow varying 
(classical) parameters set S and the mass flux of the same species J2{r, t) as a fast variable 
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on the extended set JF so that G = S [j J^. These latter variable will take into account the 
presence of inhomogeneous regions (concentration domains in this case) and so will correct 
the predictions based on the local equilibrium hypothesis. The non-equilibrium Gibbs free 
energy for a binary fluid mixture under an isotropic mechanical environment {dP = 0) then 
reads: 



The non-equilibrium contribution to the generalized entropy is given (from a Maxwell- 
like relation) [|201 by: 



By substituting equation ^ into equation ^ and using a nonequilibrium Gibbs-Duhem 
relation ll22l a dynamic contribution to the concentration of species 1 is derived namely, 




If we now assume that Antoine's law /i^(r, t) = ij) + H{r)T{f, t) holds locally, H{f) 
is a time independent parameter playing the role of the constant R (ideal gas constant) of the 
equilibrium Antoine's law; then we can derive a closed expression for for a given set of 
constitutive equations for the fast variables (in this case the mass flux and its driving force). 
We propose linear constitutive equations with exponential memory kernel for the mass flux 
J2 and its conjugated force X'^ for the following reasons: 

1. The associated transport equations are hyperbolic (of the Maxwell-Cattaneo-Vernotte 
type) Il23ll24| so causality is taken into account. 

2. These hyperbolic transport equations are compatible with the postulates of EIT lfT8ll251 



dtG = ~^dtT + iiUtC2 - ■ dth 



(4) 



Here /i^ = /i2 — j^i is the so-called relative chemical potential, J2 is the mass flux of 
species 2 and = X2 — Xi is the thermodynamic driving force conjugated to this mass 



flux. 




(5) 



|26l|27l^ 



The constitutive equations are therefore chosen to be. 




(7) 




(8) 



The Aj's are time-independent, but possibly anisotropic amplitudes, m is a unit vector in the 
direction of mass flow (the nature of u will not affect the rest of our description, since we 
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will be dealing with the magnitude of the mass flux | J2 1) and Xj's are the associated relaxation 
times considered path-independent scalars. Since we have a linear relation between thermo- 
dynamic fluxes and forces some features of the Onsager-Casimir formalism will still hold. 

By using linear constitutive equations we are limiting our treatment to comparatively 
small spatial variations, otherwise we will have to make a gradient expansion, namely J ~ 
V;U, or more generally J = Y.n ^n^^l^^ where e„ are the expansion coefficients. We believe 
that introducing the gradient expansion of the constitutive equations will not supply a big 
deal of new physics since the first contribution will be the most important from the stand- 
point of thermodynamics. We are instead considering at a deeper level the dynamic effect 
given by the exponential memory kernel. We have already talked about two fundamental ap- 
proaches in EIT, one is expanding the gradients of the non-locally conserved densities; while 
the other consists in enlarging the set of relevant parameters by adding nonequilibrium fluxes 
and forces [18J. In this work we choose to take the latter approach by considering it more 
appropriate in dealing with dynamic critical phenomena for the reasons just mentioned. We 
did just so in Eqns. [5l[6]and[8l 

As we can see, equations ([7]), ([8]) together with Antoine's law define a dynamical coupling 
between the thermal field and the concentration field. We will take the case in which the 
temperature field decays exponentially (a so-called multiplicative thermostat) since this is a 
common situation in experimental studies near the critical point 11291 . Therefore, 

T = Te, + T(f) e-"* (9) 

where T is the non-equilibrium local temperature, T^q is the temperature associated with 
the target equilibrium state, namely, the closest equilibrium state to the critical point and the 
one towards which the system relaxes in time, T(r) is an effective amplitude for the thermal 
decay rate and n is a dynamic critical exponent related to thermal relaxation whose explicit 
nature will be evidenced later by considering its relation with the distribution of decay rates 
of concentration fluctuations. Substitution of equation|9]into Antoine's law, makes it possible 
to calculate the first derivative in equation[6l namely (^), and by substituting Antoine's law 

in equation|7]and then J2 into equation[8]we are able to write down {^^)- In a similar fashion 

one can obtain (|^), etc. After substituting the derivatives and solving the integral over J2, 
the local expression for the concentration (equation (6]) reads: 



T 



dT]{ 



Hir) 



+ 



\,X2riH{r)T{r] 
Ai(l-r2ri)-l 



1 - riTi) 



-nt 



(10) 
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In order for the integrals to be convergent and physically sound, some relations arise 
between the relaxation times ri and T2 and the nonequilibrium critical exponent n, namely 
riTi < 1 and 71x2 < 1. It is through this kind of relations that an irreversible thermodynamics 
study results enlightening, since the connection between criticality and relaxation times is 
evidenced. 

3 Generalized homogeneous functions and scaling hypothe- 
ses 

3.1 Static scaling hypothesis 

In the neighborhood of a critical point one of the relevant fields (the locally conserved den- 
sities and the order parameter ^) exhibits a stronger divergence than the others. This highly 
divergent field called the order parameter field has an associated field that defines the 
generalized susceptibility. The so-called scaling hypothesis is the assumption stating that the 
dominant singularity in the free energy F(h^, T) may be rephrased in the following form (at 
least asymptotically): 



The scaling theory of static critical phenomena ^ states that the generalized susceptibility 
associated with the order parameter follows a power-law: 



is the susceptibility associated with the order parameter, ( is the correlation length 
and for the binary fluid mixture the associated universality class is that of the isotropic 3-D 
Ising system (d = 3 dimensions, n = 1 implying a scalar order parameter (concentration dif- 
ference)). For this universality class it has been found that r] = 0.04 and that C ~ (Tc — T)~^ 
with u ^ 0.62 121. 



3.2 Dynamic scaling hypothesis 

The dynamic scaling hypothesis (DSH) states that in the neighborhood of the critical point 
a scaling factor x given hy x = qC remains finite notwithstanding the divergent correlation 
length due to the fact that in the thermodynamic limit the wave number vanishes (g ^ 0). 
On the other hand, according with the theory of dynamic critical phenomena |[8l[l0l| the static 
structure factor Xq follows a power law: 




(11) 



(12) 



Xq = Te ^^(^C) 



(13) 
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with e = ^'^rjJ"^ , g{x) an unspecified monotonic function and ( the correlation length, that 
goes as C = Co^'". 

This concentration fluctuations decay rate has been related to viscous modes by Sengers 
et al lfT2l [T3l [T5l giving the dynamic scaling exponent (in Fourier space) for the non-regular 
part of the diffusion coefficient from linearized mode-coupling theory as: 

D = Doqy (14) 

where y = 1 + z^f ^ 1.0678, 2;^ being the critical exponent related to the asymptotic 
non-regular shear viscosity. 



The diffusion coefficient exhibits anisotropy D = D{q) but according with the DSH a 
scaling factor x exists such that g = | so that D = D^q^ gives D = Dq^,x remaining finite 
so may be incorporated in the amplitude Dq hence, D = DiC^^ which clearly vanishes at 
the critical point where C diverges. In this case one has: 

D = D,{Coe-'')-'^ = D^t^y = D = D2t^' (15) 

Equation (fTSl) gives us the asymptotic dependence of the diffusion coefficient with tem- 
perature in the vicinity of the critical point via the exponent /?'. By means of scaling relations 
for the non-equilibrium thermodynamic potentials (based on their property of generalized 
homogeneity) the critical exponent for this transport coefficient will be calculated. Once this 
value for [3' is known, we will be able to calculate u, since u = The value ontained for 
the exponent v could be compared with the DSH-mode coupling result in order to establish 
compatibility criteria. 



We will begin by analyzing the scaling hypothesis in the representation of the symmetric 
thermodynamic potential (/*) for a non-equilibrium binary mixture described by the time- 
dependent fields temperature (T), concentration of one species (Cj) and mass flux of this 
species (Jj): 

r = nT,Q,Ji) (16) 

Equation (fT6l) is a non-equilibrium extension of equation ([TT]) . the order parameter is 
given by the concentration field Cj, and the non-equilibrium contribution is given by the 
mass flux Jj. We will assume here that (generalized) mono-scaling applies : 

r (A'^-e, A"^C+, X-'^J^ = A/*(e, C+, Jl) (17) 

where e = = , Cc the critical consolute composition and A an arbitrary 

parameter. I and II refer to the corresponding phase. Partial derivatives of /* respect to its 
natural variables are conjugated. 
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df* ^ df* 

(^) = /.+ (e, C\ J.); (^) = -X,(e, C\ J.) 

/it - /it 

/i+ = n^ = ni-fi2. (18) 

Since /* is considered to be a generalized homogeneous function, it follows that: 

X--'-^X,{\^Ee,\''^C^,\''^Ji)=Me,c^,J) (19) 
X^^-^fi^iX^'^e, A'^^C"*", A'*-^ Ji) = /i+(e, C+, J^) (20) 
which implies the following scaling relations: 

/i+ = /io(-e)^; X = Xoi-ef 
/i+ = /io(C7+)^; X = Xi(CY' 

/i+=/io(J)^; X = X2iJY' (21) 

From now on we will obviate the subscripts in Jj and Xi and we will work with their 
scalar magnitudes since we are interested in the thermal effect of the mass flow given by 
these scalar magnitudes. Simple relations between the exponents may be derived as it is 
usual Il34l ; namely, 



(3 = ; (3 = ; 6 = ; 6 = ; r = ; r = (22) 

aE CLE 0,c 0,c 0,J CLj 

These equations may be casted to give relations between exponents related to either non- 
equilibrium or equilibrium quantities (Remember that partial derivatives as /i+ are calculated 
at constant values of J, that is steady states that can made be coincident with reversible states 
(^ S)) so we get: 



_l_r' _ _ ^ _ ^ _ ^' 

~ TT5 ~ 5'{l + r'y ~ Y+V' ~ r(l + 5)' '^'^ ~ + 5) ~ /3'(l + r') 

As we have seen equilibrium quantities are now being supplemented with a non-equilibrium 
exponent in order to determine the others, due to the fact that only 3 of these exponents are 
independent, two of them from equilibrium and one for the nonequilibrium processes. In the 
case that concerns us, we want to get relations leading to the determination of exponent (3' 
that, as is clear from Eq. [22] is given in terms of aj and aE- Eqns. [23] show that aj and 
aE are related to either the two equilibrium exponents (3 and 5 in the case of a^; or to the 
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non-equilibrium exponent r' in the case of aj. 

Under the particular approximation made in this work (namely, linear but causal pro- 
cesses) , the constitutive equation [8] implies the value r' = 1. The validity of this assumption 
in the context of dynamic critical phenomena is left for an a posteriori analysis based on the 
consequences in the final results. Nevertheless, given the dynamic character of critical phe- 
nomena, we believe that the temporal component (taken into account in the memory equation 
Eq. [22] for the nonequilibrium contribution) will be of greater importance than the spatial 
inhomogeneities. For these reasons we assume the dynamic critical exponent to have the 
value r' = 1 and relate it with the two additional equilibrium exponents in order to make the 
following calculations. 

The diffusion coefficient behaves as a power-law with the reduced temperature e with an 
exponent (3' as: 

DVc = J- J ^el^' ^ D^e^^' (23) 

Equation (l23l) has been derived from the scaling procedure of the extended thermody- 
namic potentials. If we recall equations (|22|) it has been possible to see that only three of 
the exponents are independent. Having noticed this and being aware of the values of (3 and 5 
from equilibrium thermodynamics, as well as the value of r' = 1 from our linear constitutive 
relation [equation ([8])] it is possible to find out the value of [3' and from it the value of the 
exponent v. Since we are interested in a complete thermodynamical characterization of crit- 
ical dynamic phase transitions we will also look at the behavior of the associated correlation 
function. 



4 Non-equilibrium correlation functions 

It has been proposed that macroscopic correlation functions may be obtained from a non- 
equilibrium thermodynamical description, by using the fact that extended thermodynamical 
potentials belongs to a certain class of statistical systems whose probability distribution func- 
tions obey a certain relaxation times hierarchy on its evolution towards the most probable 
state. 

We will use the equivalence between 2-time correlation functions and time evolution op- 
erators acting on a dynamical variable [35 J in order to calculate these correlation functions. 
The dynamic variables are the set of extended thermodynamical potentials. Time evolution 
operators are obtained as inner products of the aforementioned thermodynamical potentials. 
If the thermodynamic average is performed under isotropic canonical conditions [|38l[39ll40ll 
the resulting correlation function is given by: 

(Ci(f, t), Ci(f, t')) = f E{f) e-"* e"'^'*' W{t - t')dt' (24) 
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At this point it becomes necessary to clarify the meaning of the dynamic exponent n. 
Exponent n is, in principle, a time-dependent quantity characterizing the decay rate of the 
thermal field at time t and n' its value at time t'. In the approximation taken in this work 
this decay rate is closely related to the decay rate of concentration fluctuations T. For similar 
systems a distribution of decay rates r(t) has been proposed [[T21IT31 giving in some cases, a 
double exponential decay behavior for the correlation function. 

Nevertheless, in the case of a simple mixture (non-colloid, non-micellar, etc.) a simple 
exponential fits well the experimental data. This is accomplished if we consider the dynamic 
exponent n to have a weak time dependence, so that n = n' at least along the integration 
path $. In doing this we are excluding crossover phenomena, as well as some features of 
complex fluid behavior of our analysis. It is worth mentioning that recently, some complex 
fluids have been shown to follow this exponential decay of the correlation function. We 
will mention the case of an ionic binary (ethylammonium nitrate-n-octanol) mixture [41], 
and also binary mixtures of small molecular weight fluids (like hexane-nitrobenzene) [42] . 
In contrast, a protein in aqueous solution has been observed to significantly deviate form 
exponential behavior [|43l . In the single exponential approximation one has: 



(Ci(f, t), Ci(f, t')) = f S(r) e-"(*+*') W{t - t')dt' (25) 

here W {t — t') is a weighting function and the following abbreviations have been made, 

S(r) =AB + B'^ (26) 

\i\2nH{r)T{r) ^ Ain Ai(l -nri) - 1 

^= n T^Ji^i^i ~ ^2(7^ r + 7--. ^ )J (28) 

n(\ — riTiY (1 — riTi) (1 — nT2) 



Equation (125]) is the non-regular part of the composition field time correlation function. 

One of the advantages of this formalism with respect to others, such as mode-coupling 
theory is that we can test different kinds of relaxational couplings (delta-correlated, quasi- 
markovian, gaussian, lorentzian, Omstein-Uhlenbeck type, non-markovian, simultaneous 
non-linear multi-modal coupling, etc.) just by changing the weighting functions W{t — t') 
(that is the family of probability measures), within a unifying thermodynamic scheme. 



For a broad family of weighting functions, equation (|25]) asymptotically converges to a 
limit given by P4l : 

(Ci(f, t), Ci(f, t')) = K S(f) (29) 

where k is a constant depending on the explicit weighting function . For an Omstein- 
Uhlenbeck measure (W{t — t') = e'*~*'') convergence implies k = In the case of 
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Gaussian measure (W{t — t') = e*^*"*')^) an asymptotic solution similar to equation (|29l) 
could not be given in terms of an exponential decay. In this case the solution consist of 
two contributions: an exponential decay plus an error function type mode. This difference 
eliminates the short-time plateau (known as critical slowing down of fluctuations) present in 
correlation functions obtained using colored noise measures as weighting functions. 



5 Results and discussion 

Once we have been able to calculate the correlation function, we are now in position to 
discuss some interesting physical results derived from the scaling analysis of this function 
and its comparison with experimental data. 



5.1 Diffusion coefficient critical exponent 

The concentration correlation function (as obtained by light scattering or inelastic neutron 
scattering) is given by the following (van Hove) expression : 

(Ci(f,t)Ci(f,t')) = x(r)e-^'^'* (30) 

It is a known fact that the dynamics of the fluctuations slow down close to a critical point. 
If we look at the critical consolute point of a binary fluid mixture a sound physical reason 
for this critical slowing down could be extracted from a simple analysis of equation (|30l ) for 
the 2-time concentration correlation function. In this case the diffusivity D is to be identified 
with the mass diffusion coefficient for one of the species between the two phases. At the 
consolute point (where the two phases coalesce in one single phase) there is no driving force, 
for the diffusion mechanisms since spatial concentration gradients no longer exist. 

Considering an exponential thermal decay (a multiplicative thermostat) from the critical 
point [equation ©J we prepare the system at the critical point at t = with the critical 
temperature given by: 

T, = T,, + T{r} (31) 

The correlation length diverges according with the power law: 

c ~ ij^^r (32) 

In our case: 

Tc 

C = ^o(T(f)(l-e-*)^' ^^^^ 
at t=0 (critical point) T = Tc y C ^ oo. 
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In the long-time limit the system will reach an equilibrium state (the target state) with 
the following properties, at t ^ oo, T = Tg^, C = Co (5^)''-The physical meaning of this 
equilibrium state and its associated temperature and correlation length will become clearer 
later. 



Let us recall equation [25J since the critical concentration correlation function decays as 
e~^"*, for a variety of non-equilibrium measures W{t') [as stated by equation (|29l )1 we get 
the following relations: 

(Ci(f, t), Ci(f, t')) ~ ~ e-^^'* (34) 

It has been proved useful to define a decay rate of fluctuations T = Dq^ according to 
equation |30] as an indicator of the decay velocity of the concentration correlations. This 
definition clarifies the meaning of our dynamic critical exponent n since T = 2n. Following 
the behavior of the diffusion coefficient given by Sengers et al [[T3l [T4l [TSl namely equation 
(fT4)) . we have: 

g-2nt ^ g-Do,«+^t. ^ ^ £^ (35) 

but the correlation length depends on the time-dependent temperature field according to 
equation (|33l) hence, 

n=\^\{ .. I J '" (37) 



2 / \ Co( 



T(r)(l-e-"*)' 



and so. 



^ = :^^) (3^) 



Since the decay rate of concentration fluctuations V = 2n equation [38] can be rephrased 



(39) 



The critical exponent u may therefore be obtained from the thermodynamical description 
given by equation (|38l ) as follows: 

zy=lini(^— )^ L (40) 
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Equation |40] is the main result of this paper, by means of it and by multiplying by 
y = 1 + z,, a value for P' is obtained to compare with the value of equation [151 The value for 
u was obtained by numerical evaluation of the limit (lim^^ q). We made this by evaluating the 
function for several values of t and then by imposing an error bound of 0.001 in the quantity 
calculated. We looked up for the accumulation point of the function near its zero-time limit, 
that is, we use a Cauchy convergence criteria in the limit. This leads us to adjust a cut-off time 
value of the order of 95 milliseconds. Our closest approach to the limit given by this value 
of t. Given this facts, we should expect an error bounded in the third significant decimals of 
exponent u since Cauchy criteria is linear on the error bound. 

It is possible to rearrange the preceding expression (Eq. HO]) to give: 



{x/CoY^y+^\i-e"''Y^y+^^ (4i) 



We now define fi, and call it the modified dynamic Kawasaki function for reasons to be 
seen, as follows: 

Substitution of equation |4T] into equation |42] yields finally that, 

m = (a^/Co) (1 - e-^y (43) 



We wrote as a function of time via the decay rate n and not just as a function of the 
scaling parameter x in order to avoid the use of the Ornstein-Zernike (OZ) approximation. 
We notice that Vt is finite and regular at t = (Critical point) and in the t oo limit 
(thermodynamic equilibrium). 

n{t = 0) = 0; n{t ^ oo) = (x/Co) (44) 

This modified Kawasaki function (scaled by the exponent Lp = possess a clearer 
physical interpretation than the one obtained from the asymptotic Ornstein-Zernike form of 
the radial distribution function g(x) = as treated by means of Mode Coupling Theory 
for the case of non-equilibrium states since the closure of the OZ equation is given by a local 
equilibrium distribution. 



If we take experimental values for the amplitudes (Dq, T{f)) and other characteristic pa- 
rameters [42, 45], it is possible to express a convergent sequence leading to the determination 
of u from equation ([401) . In order to calculate u we had to adjust a cut-off time of the order 
of 95 ms (or a cut-off frequency of 10.52 Hz). To within this approximation we obtain an 
average value of Ucaic = 0.6303 for the set of realizations for the limit in equation[40l Once 
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this critical exponent is known we proceed to express the behavior of other physically rele- 
vant quantities. 

In fig. 1 we plotted the predicted behavior of the mass diffusion coefficient D in the prox- 
imity of the critical point. According with the DSH treatment in the already mentioned mode- 
coupling description [fT2| D = D2 e^^ = D2e^ , if we use our calculated value Ucaic = 0.6303 
we get the solid line plotted in fig. 1. We also plotted experimental results for an ionic binary 
mixture as reported from recent scattering experiments [|45l . As is shown good agreement be- 
tween experiment and theoretical predictions is achieved even for this case (an ionic mixture). 
We must anyhow say a word of caution about the effect of using a reduced temperature coor- 
dinate since this fact stretches the curve and thus gives a visual impression of higher accuracy. 
What is surprising is that even for this relatively simple model (that leaves out many char- 
acteristic features) essential behavior is reproduced: a signature of a thermodynamic analysis. 

In fig. 2 the wave number dependency of the diffusion coefficient is shown. From some 
time ago experimentalists have claimed a non-fickian behavior near the critical point against 
some theoretical assumption of linearity derived from generalized hydrodynamics. Fig. 2 
shows that even though this deviation from Fick's law exists indeed the effect is really small. 
Our formalism predicts D ~ with y = 1.0678. Since this was obtained using the mono- 
scaled DSH it is also possible that further corrections need to be taken into account in the 
future. This fact will not affect significantly our consideration made in equations [8]regarding 
r' = 1. 

In fig. 3 we show that the decay rates of concentration fluctuations deviates slightly from 
the parabolic behavior given by linear response predictions. Fig. 3 shows a semi-flat region 
in the low wave number regime that possibly indicates a macroscopic effect of stability at 
long distance scales. 

In Table 1 we compare our results with some of the most reliable data published in the 
recent literature of dynamic critical phenomena. The value of exponent u is somewhat higher 
than the reported values. This fact may be explained because the experimental data we used 
for our calculations were measured on a region that, appart from critical phenomena, ex- 
hibits other complex dynamic processes like Rayleigh-Benard convection due to temperature 
fluctuations. Our calculation of u is based on the behavior of the decay rate of thermal fluc- 
tuations r according to equations [39] and |40l A cross-over in the associated structure factor's 
dependence of the wave number q is known to exist on passing from large scale to small scale, 
and this fact has been related with the presence of two modes in the Rayleigh spectrum, one 
thermal mode identified with thermal fluctuations of decay rate DtQ^ and one viscous mode 
with a decay rate Ag^ [|5TI . Here Dt is the thermal diffusivity and A is the kinematic viscosity. 
In our theoretical approach we did not take viscous effects into account except in the con- 
sideration for exponent y = 1 + 2;^ in equation [14] and so viscous dissipation is disregarded 
in our treatment of the dynamics of thermal fluctuations. Since the viscous mode is faster. 
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the difference in the dynamics between this calculation and others should vanish eventually, 
although an effect in the critical exponent, even if small, is present. 

An additional analysis of fig. 3 suggest that it is possible that, just as there is a dynamic 
effect of slowing down of critical fluctuations, there is also a spatial effect of smoothing of 
critical fluctuations. Further investigations are need in the future in this direction and, possi- 
bly by inclusion of other dissipative effects (viscous modes notably). Our physical intuition 
about the applicability of non-equilibrium thermodynamic description to critical behavior 
will be substantially increased. 
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Figure captions 

Figure 1 Temperature dependence of the near critical mass diffusion coefficient. Solid 
line is this works prediction [equation (|23l) 1. Square dots and open circles are experimental 
values taken from reference [|45l . 

Figure 2 Wave number dependence of the near critical mass diffusion coefficient. A 
slight deviation of Fickt's law could be noted. The curve represents a power-law behavior 

D ~ gl 0678_ 

Figure 3 Decay rate of fluctuations F as a function of the wave number. This work 
predicts a slight deviation from the linear response parabolic behavior. It is noticeable the 
presence of a plateau in the low wave number region which possibly represents a macroscopic 
effect of spatial smoothing down of fluctuations. 
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Critical exponent 


Theoretical calculations 


Renormalization & simulations 
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0.6303 ± 0.001 This work 




V 


0.63012(16) HTS [461 


0.618 |47| 


V 


0.6302(2) HTS [461 


0.630 ± 0.00015 |49| 


V 


0.63002(23) HTS [47 1 


0.632 ±0.003 |49| 



Como vemos en la tabla ?? esto es chido. 
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Figure 1: Temperature dependence of the near critical mass diffusion coefficient. Solid line is 
this work's prediction [equation (l23l)1. Square dots and open circles are experimental values 
taken from reference ll45l 
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Figure 2: Wave number dependence of the near critical mass diffusion coefficient. A slight 
deviation of Fickt's law could be noted. The curve represents a power-law behavior D ~ 

„ 1.0678 
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Figure 3: Decay rate of fluctuations F as a function of the wave number. This work predicts 
a slight deviation from the linear response parabolic behavior. It is noticeable the presence of 
a plateau in the low wave number region which possibly represents a macroscopic effect of 
spatial smoothing down of fluctuations. 




